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, i ' Abstract. The pseudo-rapidity dependence of the elliptic flow at various excitation energies mea- 

sured by the PHOBOS Collaboration in Au+Au collisions at RHIC is one of the surprising results 
that has not been explained before in terms of hydrodynamical models. Here we show that these data 
are in agreement with theoretical predictions and satisfy the universal scaling relation predicted by 
the Buda-Lund hydrodynamical model, based on exact solutions of perfect fluid hydrodynamics. 
We also show a theoretical prediction on the rapidity and transverse momentum scaling of the HBT 
\ radii measured in heavy ion collisions, based on the Buda-Lund model. 
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INTRODUCTION 

O 

One of the unexpected results from experiments at the Relativistic Heavy Ion Collider 
(RHIC) is the relatively strong second harmonic moment of the transverse momentum 
distribution, referred to as the elliptic flow. Measurements of the elliptic flow by the 
PHENIX, PHOBOS and STAR collaborations (see refs. [1, 2, 3, 4]) reveal rich details 
in terms of its dependence on particle type, transverse and longitudinal momentum 
variables, on the centrality and the bombarding energy of the collision. In the soft 
transverse momentum region, these measurements at mid-rapidity are reasonably well 
described by hydrodynamical models [5,6]. However, the dependence of the elliptic flow 
on the longitudinal momentum variable pseudo-rapidity and its excitation function has 
resisted descriptions in terms of hydrodynamical models (but see their new description 
by the SPHERIO model [7]). 

Here we show that these data are consistent with the theoretical and analytic predic- 
tions that are based on eqs. (1-6) of ref. [8], that is, on perfect fluid hydrodynamics. 

We furthermore calculate rapidity dependent HBT (Bose-Einstein) radii in the frame- 
work of the model and make prediction on the universal scaling of these observables. 

Our tool in describing the pseudorapidity-dependent elliptic flow and HBT radii is 
the Buda-Lund hydrodynamical model. The Buda-Lund hydro model [9] is successful 
in describing the BRAHMS, PHENIX, PHOBOS and STAR data on identified single 
particle spectra and the transverse mass dependent Bose-Einstein or HBT radii as well 
as the pseudorapidity distribution of charged particles in Au + Au collisions both at 
-y/sNN = 130 GeV [10] and at ^/snn = 200 GeV [11]. However the elliptic flow would 
be zero in an axially symmetric case, so we developed the ellipsoidal generalization of 
the model that describes an expanding ellipsoid with principal axes X, Y and Z. Their 
derivatives with respect to proper-time (expansion rates) are denoted by X, Y and Z. 



The generalization goes back to the original one, if the transverse directed principal 
axes of the ellipsoid are equal, ie X = Y (and also X = Y). 

The deviation from axial symmetry can be measured by the momentum-space eccen- 
tricity, 

The exact analytic solutions of hydrodynamics (see ref. [8, 12, 13]), which form the 
basis of the Buda-Lund hydro model, develop Hubble-flow for late times, ie X — Xx, 
so the momentum- space eccentricity e p nearly equals space-time eccentricity £. 

Let us introduce At] additionally. It represents the elongation of the source expressed 
in units of space-time rapidity. Let us consider furthermore that at the freeze-out xAr] = 
Z and Z Zx, and so Ar\ m Z 

Hence, in this paper we extract space-time eccentricity (£), average transverse flow 
(u t ) and longitudinal elongation (Ar|) from the data, instead of X, Y and Z. 

In the time dependent hydrodynamical solutions, these values evolve in time, however, 
it was show in ref. [14] that X, Y and Z, and so £, u t and Ar] become constants of the 
motion in the late stages of the expansion. 



RAPIDITY DEPENDENT ELLIPTIC FLOW 

The result for the elliptic flow (under certain conditions detailed in ref [15]) is the 
following simple universal scaling law: 

Iq{w) 

The model predicts an universal scaling: every V2 measurement is predicted to fall on 
the same universal scaling curve I\/Iq when plotted against w. 

This means, that V2 depends on any physical parameter (transverse or longitudinal 
momentum, center of mass energy, centrality, type of the colliding nucleus etc.) only 
through the (universal) scaling paremeter w. 

Here w is the scaling variable, defined by 

W = &(±-f), (3, 

Am t \T* ; y T* !X J 



and 



T*, x = T + m t X 2 T °_ (4) 
lo + m t a z 

T^ y = T + m t Y 2 (5) 



and 



m t =m f cosh(r|i-y). 



(6) 



Here a = (AT /T) t measures the temperature gradient in the transverse direction, at the 
freeze-out, m t is the transverse mass, To the central temperature at the freeze-out, while 
r\ s is the space-time rapidity of the saddle-point (point of maximal emittivity). This 
saddlepoint depends on the rapidity, the longitudinal expansion, the transverse mass and 
on the central freeze-out temperature: 



r\ s -y 



y 



1+Arj^ 



(7) 



To 



where y = 0.5 log is the rapidity. 

More details about the ellipsoidally symmetric model and its result on V2(t|) can be 
found in ref. [15]. 

Eq. 2 depends, for a given centrality class, on rapidity y and transverse mass m t . 
Before comparing our result to the V2(t|) data of PHOBOS, we thus performed a saddle 
point integration in the transverse momentum variable and performed a change of 



variables to the pseudo-rapidity T] = 0.51og( 



\P\+Pz - 
\P\-Pz- 



similarly to ref. [16]. This way, 



we have evaluated the single-particle invariant spectra in terms of the variables r\ and (]), 
and calculated V2(t|) fr° m this distribution, a procedure corresponding to the PHOBOS 
measurement decribed in ref. [1]. 

We have found that the essential fit parameters are £ and At], and the quality of the fit is 
insensitive to the precise value of To, a and u t . These parameters dominate the azimuthal- 
averaged single particle spectra as well as the HBT (Bose-Einstein) radii, however they 
only marginally influence v%. Their precise value is irrelevant in a broad region of values 
and does not influence the confidence level of the V2(t|) fits. Hence we have fixed their 
values as given in the caption of table 1 . We also excluded points with large rapidity 
from the fits in case of lower center of mass energies. 

Fits to PHOBOS data of ref. [1] and its l-3o error contours are shown on the top 
two panels of fig. 1. The fitting package is available at ref. [17]. Bottom panel of 
fig. 1 demonstrates that the investigated PHOBOS data points follow the theoretically 
predicted scaling law. 

Table 1. Results of fits to PHOBOS data of ref. [1]. Both space-time eccentricity 
(e) and longitudinal elongation (Ar|) increase with increasing y'SNN- Remaining 
parameters were fixed as follows: Tq = 175 MeV, a = 1.19 and u t = 1.64. 





19.6 GeV 


62.4 GeV 


130 GeV 


200 GeV 


e 


0.294 ± 0.029 


0.349 ± 0.008 


0.376 ± 0.005 


0.394 ±0.006 


AT) 


1.70 ±0.25 


2.16 ±0.05 


2.46 ± 0.04 


2.56 ± 0.04 


X 2 /NDF 


1.84/11 


20.1/13 


34.8/15 


27.5/15 


conf. level 


100% 


21.4% 


1.00% 


7.03% 




Figure 1. Top: PHOBOS data on the pseudorapidity dependence of the elliptic flow [1], at various 
center of mass energies, with Buda-Lund fits. Middle: Error contours of the fits. Bottom: Elliptic flow 
versus scaling variable w is plotted. The data points of ref. [1] show theoretically predicted [15]) universal 
scaling, when plotted against the universal scaling variable w. 
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Figure 2. Upper panel: R out , R s ^ e and R\ on g as a function of rapidity y and transverse mass m t . Lower 
panel: The two-dimensional R(m t ,y) functions are predicted to show a scaling behavior, insofar as they 
depend only on scaling variable W t . 



RAPIDITY DEPENDENT HBT RADII 

In the framework of the model we can also calculate the HBT radii. In the simplest 
case, where system of ellipsoidal expansion equals the out-side-longitudinal coordinate 
system: 

*U - z^i+^lE^V'. do) 

This means, that the HBT radii depend on transverse mass and rapidity only through 
the scaling paremeter m~ t , as illustrated on fig. 2. This behavior could easily be checked 
by measurement of rapidity and transverse momentum dependence of the HBT radii 
and comparing this data to the present prediction of the Buda-Lund model. Such a 
comparision could be a further test of perfect fluid hydrodynamics. 



CONCLUSIONS 



In summary, we have shown that the excitation function of the pseudorapidity depen- 
dence of the elliptic flow in Au+Au collisions is well described with the formulas that 
are predicted by the Buda-Lund type of hydrodynamical calculations. 

We have provided a quantitative evidence of the validity of the perfect fluid picture of 
soft particle production in Au+Au collisions at RHIC but also show here that this perfect 
fluid extends far away from mid-rapidity. 

We also suggest a further test of perfect fluid hydrodynamics at large rapidities, 
expressed by eqs. 8-10 and illustrated by fig. 2. 

The universal scaling of PHOBOS V2(t|), expressed by eq. 2 and illustrated by fig. 1 
provides a successful quantitative as well as qualitative test for the appearence of a 
perfect fluid in Au+Au collisions at various colliding energies at RHIC. 
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